Deconfined SU(2) phase with a massive vector boson triplet 
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We introduce a model of SU(2) and U(l) vector fields with a local U(2) symmetry. Its action 
can be obtained in the London limit of a gauge invariant regularization involving two scalar fields. 
Evidence from lattice simulations of the model supports a (zero temperature) SU(2) deconfining 
phase transition through breaking of the SU(2) center symmetry, and a massive vector boson triplet 
is found in the deconfined phase. 
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I. INTRODUCTION 

Euclidean field theory notation is used throughout this 
paper. The action of the electroweak gauge part of the 
standard model reads 



S 



L cw = 
p™ = 
F b = 



d A xL cw , (1) 

_1 pompom rry T?b pb (2) 

d^CLu - dvdfj, , (3) 

d^B v -d v B^ + ig b [B^B u ] , (4) 



where a'^ are U(l) and B^ are SU(2) gauge fields. 

Typical textbook introductions of the standard model, 
e.g. [l|, emphasize at this point that the theory con- 
tains four massless gauge bosons and introduce the Higgs 
mechanism so that one obtains a massive vector bo- 
son triplet and only one gauge boson, the photon, stays 
massless. Such presentations reflect that the introduc- 
tion of the Higgs particle in electroweak interactions 
|H preceded our non-perturbative understanding of non- 
Abelian gauge theories. 

In fact, massless gluons are not in the physical spec- 
trum of (JT]). The self interaction due to the commutator 
(U]) generates dynamically a non-perturbative mass gap, 
and the SU(2) spectrum consists of massive glueballs. 
The lightest glueball can be used to set a mass scale. 
Choosing for it, e.g., 80 GeV and coupling fermions is 
perfectly admissible. This does not constitute an ansatz 
for an electroweak theory by two reasons: The SU(2) 
glueball spectrum is not what is wanted (e.g., masses of 
spin and 2 states are lower than for spin 1 Q) and 
fermions would be confined into boundstates, which is 
not the case. Coupling a Higgs field causes a deconfining 
phase transition, so that fermions are liberated, a pho- 
ton stays massless and glueballs break up into elementary 
massive vector bosons. Such a conhnement-Higgs transi- 
tion has indeed been observed in pioneering lattice gauge 
theory (LGT) investigations 4]. The present paper in- 
troduces a different model for which a zero-temperature 
deconfined phase breaks the Z2 center symmetry of SU(2) 
and exhibits a massive vector boson triplet. 

The motivation for the action defined in the next sec- 
tion comes from properties of the U(l) Polyakov loop. 
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FIG. 1: Scatter plot for U(l) Polyakov loops on a 12 4 lattice 
at p a = 0.9 in the symmetric phase (center) and at /3 a = 1.1 
in the broken phase (ring). 



U(l) LGT confines fermions in the strong coupling limit 
of its lattice regularization [5| . At weaker coupling it un- 
dergoes a transition into the Coulomb phase as has been 
rigorously proven (f|. In the Coulomb phase the effec- 
tive potential of the U(l) Polyakov loop P a assumes the 
Mexican hat shape that is also characteristic for a com- 
plex Higgs field in the broken phase. This is illustrated in 
Fig. [T]by scatter plots of Polyakov loops in the symmetric 
and in the broken phase. From the ring-like distribution 
in the broken phase the similarity with the behavior of a 
Higgs field is evident (details of the simulations for the 
figure are given in section! 



We introduce a contribution to the action that is the 
alternating product of SU(2) and U(l) gauge matrices 
around a plaquette. It unifies SU(2) and U(l) gauge 
transformations into a local U(2) symmetry, which we 
take as the starting point of our presentation in sec- 
tion mi a gauge invariant regularization is obtained by 
coupling two scalar fields and leads to our model by 
freezing these fields up to their gauge transformations, 
here called London [§| limit. The scalar fields can then 
be absorbed by extended gauge transformations, which 
were introduced and discussed in 0, [l(| ■ Some techni- 
cal details of section[H]encountered in the investigation of 
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the classical continuum limit are relegated to appendixes. 

In section IIIII we investigate the quantum field theory 
numerically in the lattice regularization. Evidence for a 
SU(2) deconfming phase transition, first observed in 
and results for the particle spectrum are presented. A 
massless photon is found in both phases. Noisy correla- 
tions (as usual) indicate massive SU(2) glueballs, which 
appear to jump to unphysically high mass values in the 
deconfined phase. A beautiful strong signal for correla- 
tion functions of a massive vector boson triplet is seen 
in the deconfined phase, while it does not exist on the 
confined side. 

Summary, outlook and conclusions follow in the final 
section |IV] 

II. ACTION AND GAUGE 
TRANSFORMATIONS 

Introducing gauge fields is often motivated by promot- 
ing a global symmetry of matter fields to a local one. Let 
us assume a global U(2) symmetry: 

ip -> ip' = Gip with G = GiG 2 , (5) 
Gi G U(T) and G 2 G SU{2) . 

The ipj^d^tp part of the free fermion Lagrangian breaks 
([5]) as a local symmetry. This is overcome by coupling to 
gauge fields, whose transformation behavior cancels the 
unwanted contributions out. In the lattice regularization 
the relevant coupling reads 

a 3 K fi ip(na) J7 M (na)V^ l (na)'!/;(na + fia) , (6) 

where n — (m, n^) is a vector of integers, fi is the 

unit vector in \i direction, a the lattice spacing (x — na), 

Uf, = e i3aaA " G 17(1) and = e igbaB " G SU(2) . (7) 

Here g a and gb are introduced as U(l) and SU(2) coupling 
constants. The precise definition of the fermions (see 
for instance [n], [l2j|) in Eq. ([5]) is accomplished by 
and irrelevant in our context as we discuss only their 
gauge transformations. The classical continuum limit is 
for a — > obtained from a 4 contributions, whose sum 
becomes an integration over space. 

In the lattice formulation it is obvious that the local 
gauge transformations 

ip(na) — > G{na) tp{na) = G\{na) G2{na) ip{na) (8) 

will be absorbed by imposing the transformations 

U^(na) — > Gi(na)U IJ- (na)Gi 1 (na + fia) , (9) 
Vfi(na) -> G 2 {na)Vi 1 (na)G 2 ~ 1 (na + fia) , (10) 

on the vector fields, provided that their contributions to 
the action are gauge invariant. This is (for instance) 
fulfilled for the Wilson action 

S sauge = ^^ReTrC/ p + ^^Try p (11) 
p p 



where the sums are over all plaquettes and U p , V p are 
oriented products of gauge matrices around the plaquette 
loop. E.g., for a V p plaquette in the \w, fi ^ v plane: 

V^(x) = Tr [V M {x) V„{x + fia) V}(x + va) Vj(x)] . 

As it turns out to be convenient, we represent here and 
in the following U(l) phase factors by diagonal 2x2 
matrices with a^ (JT]) and ([7]) related by 

^ = ^ofl,, (12) 

where To is the 2x2 unit matrix. The electromagnetic 
contribution to the Lagrangian ([1]) reads then 

L a = -\^F^F^, = dpAv-d^. (13) 
The equations 

1 4 

Pa = — and b = -2- (14) 

9a 9b 

relate (3 a and @b of (fTTj) to the couplings g a and gb of . 
More lattice actions that give the classical continuum 
limit ([T]) are discussed in jl3l |. 

There is an asymmetry in the way the gauge transfor- 
mations work: A transformation of either the U(l) or the 
SU(2) matrices implies the corresponding transformation 
for the fermion field, while there is no such effect on the 
other gauge field. Another remarkable [lH feature of the 
action (|11[) is that it is additive and not multiplicative 
in the U(l) and SU(2) fields. Changing also the fermion 
terms ([6]) to additive 

a 3 K fl ip(na)U /1 (na)ip(na + fia) (15) 
+ a^K^ ip{na) V fJl (na)tp(na + fia) , (16) 

the a 4 contributions to the classical a — > limit re- 
main unchanged, while the gauge transformations be- 
come symmetrical in all fields, because Eq. ([9]) and (fT0|) 
have to be replaced by local U(2) transformations 

U^na) G(na)U l +(na)G~ 1 (na + fia) , (17) 
V^(na) -> G(na)V fl (na)G^ 1 (na + fia) , (18) 

which we call extended gauge transformations while re- 
serving the notation gauge transformations for their con- 
ventional version. 

The gauge part (llip of the action is invariant under 
local U(2) transformations. After a generic transforma- 
tion former U(l) and SU(2) matrices will both be in U(2) 
and their identification as U(l) or SU(2) has become con- 
verted into constraints on the gauge manifold. The spe- 
cial gauge in which they are actually in U(l) and SU(2), 
i.e. Eq. ([TJ holds, will be called diagonal gauge. In this 
gauge the U(l) part is not only diagonal but proportional 
to the unit matrix. With reference to the existence of this 
gauge we keep the U(l) and SU(2) matrix identifications. 
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Requiring invariance only under extended gauge trans- 
formations p7|) and (fT5|) allows for new additions to the 
action. Numerically we investigate 



cadd \ r*add oadd 



(19) 



-ReTr [U»(x)V v (x + jla) Ul(x + Pa) V}(x)] . (20) 
A gauge invariant regularization is [3] 

ss = Y,\J2 s U+ kT * [(^-^m, (2i) 

S% = ^ReTr {U^x)Ul(x + va) (22) 

[<&t(a; + (id) V v (x + (id) $>(x + jla + va)] 
[&(x)Vu(x)<f>(x + £>a)}^} 

where $ is a 2 x 2 matrix scalar field that is charged with 
respect to J7(l) and SU(2). Its gauge transformations 
are 



g§. 



(23) 



where 5 £ SU(2), e JQ € U(l). In the London |8j limit 
k — y 00 this action is equivalent to (fH?|) . because $ takes 
on it vacuum value, which is a pure gauge 



.9- 



(24) 



Fixing the gauge to $ = ro yields (jl9|) and on finite lat- 
tices properties at sufficiently large n values are expected 
to be practically identical. 

It is only in the limit k — ¥ 00 that the coupling con- 
stant A s together with the dimensions of the scalar fields 
becomes the dimensionless coupling A of (| 19[) (recall that 
the dimensionless lattice field $ is related to the physical 
scalar field $ p h ys by $ = a3> P hys [HI)- This is opposite 
to what happens in the SU(2) Higgs model where the 
interaction between scalar field and SU(2) field, as for 
instance given in [TBj . is proportional to 



X fJ,— l 



(25) 



with <p — pg, p > 0, g € SU(2). Here A^ is dimensionless, 
and fixing the gauge in the London limit to <j> = to leads 
to the gauge symmetry breaking contribution 



^EE TrV ^ 

X jU=l 



(26) 



where Aft acquires the dimension of the </>'</> scalar fields. 

The extended gauge transformations (fT7|) and (fT5|) re- 
store the transformation (|2"3")l without using scalar fields 
by absorbing the exp(— ia) phase factor into the and 



the SU(2) transformation g into the field. This re- 
quires the action (|19[) . or similar, and is not possible in 
the London limit (f26|) of the Higgs model. 

Before investigating the quantum field theory (|19|) nu- 
merically, we consider the classical continuum limits of 
(EGO) and (ED- 



A. Classical continuum limit 

The classical Lagrangian is expected to be the effective 
theory at energy scales much larger than the masses of 
the model. In the limit a — > the a 4 contributions of 
(EEO) give 



t add ^ rri / TTiadd TrtaddN 



(27) 



^ d = gbd^B,- g a d u A^+ig a g b [A^,B u } , (28) 

where the commutator reflects that and A„ will in 
general not commutate unless we are in the diagonal 
gauge, defined in-between Eq. (fT8|) and (CEO). 

Let us translate the extended gauge transformations 
(fTT|) and (|T5|) into continuum notation. Consider the 
gauge-covariant derivative 



(29) 



of an electromagnetic field a^(x) on a complex fermion 
field ip(x). With the gauge transformations 



1 

a>n(x) -> a' = a M 9 p a(x) 



9a 



one finds 

so that the Lagrangian 



(30) 
(31) 

(32) 



i c 7 = ff, (i llx D™ -m)l>- jF^F™ (33) 



is gauge invariant, where 



TTicm ^ [ /-) em r-jeml 



(34) 



is the field tensor. 

Assume now that ip{x) is a complex doublet, which 
transforms under local U(2) transformations 



ij -> V/ = G(a;)V, 
G(x) = exp I -E 7 "'"*^) 

V i=0 / 



(35) 
(36) 



Here To is as in (fl"2"|) the 2x2 unit matrix and tj, i = 
1, 2, 3 are the Pauli matrices. We still can couple ?/> in an 



4 



invariant way to an electromagnetic field. We define the 
gauge covariant derivative by 



(37) 



and the extended gauge transformations of by (com- 
pare for instance Ref. [l[ for usual SU(2) gauge transfor- 
mations) 

A^ —> A'^ — GA^G- 1 + — (dpG) G- 1 (38) 

9a 

which adds to the U(l) field the transformation of a null 
SU(2) held and yields the desired result 



The field tensor defined by 

f 



transforms as 



pa jp a Qpa q- 



(39) 



(40) 



(41) 



so that the Lagrangian 



so that L add is invariant. The algebra for (|46|) is given in 
appendix [A] 

As it is tedious and not very enlightening to calcu- 
late the classical continuum limit of the action (|2"Tj) with 
scalar fields, we relegate its discussion to appendix [B] 
The essence is caught by a much simpler gauge invariant 
classical interaction with two scalar fields, which reduces 
also to (|!2T|) in the London limit, but has a lattice regu- 
larization that differs from (|2lj) . 

Let us introduce the scalar fields 



p e la To, p > and <fi 2 — a g, a > 



(47) 



with g E SU(2). It is then easy algebra to show that the 
field tensor 



S 



fMU 



4DlDU 2 + (Dl^Dlfo 



(48) 



reduces by gauge fixing in the London limit to i F^ dd 
in the diagonal gauge, so that 



L s = --Tr 
4 



it 



Sf. 



becomes a gauge invariant extension of (|27[) . 



(49) 



L% = $ (i7 M £>2 - m) $ - -Tr (F^Fp 



(42) 



stays invariant. The local U(2) transformations do not 
destroy the fact that A^ describes just an electromagnetic 
field. Any A^(x) field is gauge equivalent to one in the 
diagonal gauge for which A^ is proportional to the unit 

matrix. Consequently, F^ stays always diagonal and the 
G matrices can be omitted in (|4"Tj) . 

Similarly, we extend gauge transformations of a SU(2) 
field B^(x) by a phase. The gauge covariant derivative is 

D h l = d„ + ig b B ll with B M = ~f-6 M (43) 
and the extended gauge transformations of are 

B„ -> fl' = GB^G- 1 + — (d^G) G" 1 . (44) 
9b 



Equations (|39[) to (1421) carry simply over by replacing all 
labels a by b. 

An electroweak Lagrangian of the type 



r ab 



--Tv(F a F a ) - -Tr (F b F b ) 



" ^ (45) 
+ 1> {il^F>l ~m)tP+^ (i^Dl - m) ip 



allows one to add ([27)) . Under extended gauge transfor- 
mations {38]) for A^ and (|44j) for B^ the F* dd tensor 
([28| transforms according to 



padd 



F, 



'add 



fJLU 



GTTiadd f~i — 1 



(46) 



III. NUMERICAL INVESTIGATION 

Pure U(l) LGT with the Wilson action has at j3 a 
1.01 a phase transition from its confined into its Coulomb 
phase. Presumably the transition is weakly first order as 
first reported in [lj|. We have calculated lattice aver- 
ages of U(l) Polyakov loops P a by Monte Carlo (MC) 
simulations on a 12 4 lattice with periodic boundary con- 
ditions with a statistics of 10,000 sweeps for equilibration 
and 160,000 sweeps with measurements. Measurements 
are plotted every 20 sweeps. Figure [T] compares scatter 
plots at p a = 0.9 in the disordered confined phase and at 
(3 a = 1.1 in the ordered Coulomb phase. In the unbro- 
ken disordered phase the values scatter about zero, while 
they form a ring in the broken ordered phase. As the 
transition happens at zero temperature, this holds on a 
symmetric lattice for Polyakov loops winding in any one 
of the four directions through the torus (ordered starts 
are used to avoid metastabilities of the MC algorithm). 

In the following MC simulations with the plaquette 
action 



s = S^ usc + S 



'add 



(50) 



are performed, where S ,gaugo is defined by JTTJ and S' add 
by (fT9l . We use lattice units a = 1 in this section. 



A. Integration measure and Monte Carlo updating 

Our MC procedure proposes the usual U(l) and SU(2) 
changes. For the update of a U(l) matrix U^n) we need 
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the contribution to ([50]) , which comes from the eight sta- 
ples containing this matrix 

+ Ut(n + fi-D)Ul(n-v)U v (n-v)] 

+ ^J2\ v ^ n+ ^ u t( n + i> ) v ^ a ) 

V 

+ v2{n + [i-i>)Ul{n-v)Vv{n-P)) (51) 
and correspondingly for the SU(2) matrix V^in) 

Vu,M = yJ][K(r. + A)^(n + i))Vj(n) 
+ Vj (n + /t — i>) (n - fr) K (n - *>)] 

+ 

+ Ut{n + fi-u)Vl{n-y)U v {n-u)} . (52) 

Updates are then performed with the Gibbs-Boltzmann 
weights 

dU^n) exp { ReTr [U^n) U u A n )] I . (53) 
dV„ (n) exp { ReTr [V^ (n) Kj,„ (n)] } . (54) 

where the dU^{n) integration is from — 7r to +7r over 
the phase 4>^{ n ) of the matrix U u {n) = exp[i0 Al (n)] and 
dV^n) is over the SU(2) Hurwitz [17| measure. This is 
well suited and done with the biased Metropolis-heatbath 
algorithm [l8| . 

It is instructive to discuss (|54"|) for the case in which 
the U(l) matrices in (|52")) are aligned as it is approxi- 
mately the case in the U(l) Coulomb phase. Then the 
U(l) matrices cancel out, so that large values of A fa- 
vor Tr[V^(n)Fj(n + v)] = Tr[V^(n)V^(n - 0)] = 2 and, 
hence, the SU(2) matrices are also aligned, V^(n) = 
Vfj,(n + 0) = Vn(n — 0). SU(2) deconfinement (break- 
ing of the Z2 center group) is achieved when this effect 
(spoiled by U(l) fluctuations at finite A) becomes strong 
enough. Integrations with the measures (I5U1) and ([54]) 
do not include extended gauge transformations (1171) and 
(fT8|). Updates stay within the diagonal gauge defined in 
section HU Within the MC procedure only global SU(2) 
transformations V^(n) —> GV li (n)G^ 1 with the same G 
for all SU(2) matrices remain. 

The partition function of the MC calculation 

r 4 

Z{G(n)} = / 1] II dC/ » dV ^ n ) e?W*V>' v *Vn 

(55) 

is invariant under extended U (2) gauge transformation 
(fT7]) and (fT8)) . The Jacobian determinants of both 

U.nH l i dU ^( n ) and ELII/i^M are one and tne ac- 
tion ((511)) is by construction invariant. U(l) and SU(2) 



Wilson loops are invariant operators as the U(2) trans- 
formations drop out in the trace. As usual [HI El] their 
physical interpretation can be derived by imagining a 
coupling to static quarks. In addition to the normal Wil- 
son loops traces of corresponding mixed products of U(l) 
and SU(2) matrices are also invariant, but we make no 
use of them in the following. 

With normalization of the integration over the U(2) 
measure to one, the identity 

Z = I 'jJdG(n) Z{G{n)} = Z{G{n)} (56) 

J n 

holds under extended gauge transformations. This inte- 
gration can easily be added to the updates: At site n 
all matrices on links emerging at n will be transformed 
according to 

U„ (n) ->■ G(n) U M (n) , (n) ->■ G(n) V» (n) (57) 

and all matrices on links ending at n according to 

U^n-ft) -> U^n- p)G~ l {n) , 
V^n-fi) ->■ Vp(n-pL)G- x {n). (58) 

The acceptance rate of such updates is 100% as no change 
in the value of the action is implied. 

The unusual feature is that the integration measure 
of the vector fields U^(n) and V v (n) in the functional 
integral (|55[) is not identical with the gauge measures 
in (|56p . However, for gauge transformations of scalar 
and fermion fields we are used to this and there appears 
no strong argument against having a distinction of these 
measures also for vector fields (calling them gauge fields 
may have added to the expectation that functional and 
gauge measures are the same). 

B. Zero temperature SU(2) deconfining transition 

We keep the U(l) coupling at (3 a — 1-1 and for the 
SU(2) coupling the values /3& = 2.3 is used. At A = 0, 
without interaction, (3 a is in the U(l) Coulomb phase and 
Pb in the SU(2) scaling region. Runs are performed on 
12 4 lattices. Data points are based on a statistics of at 
least 2 10 sweeps without measurements and subsequently 
32 x 2 10 sweeps with measurements. Error bars are cal- 
culated with respect to the 32 blocks. 

With increasing A one finds a strong first order phase 
transition, which is illustrated in Fig. [5] for the expec- 
tation values of the various plaquette actions. From up 
to down: (ReTrt/ p )/2 for U(l), (TrU p )/2 for SU(2) both 
contributing to QI) and X' 1 Y,^( s ^ d ) / 16 for Mixed 
(|19|). which drives the transition. Disordered starts are 
marked by "d" and ordered starts by "o" . The values are 
compiled in table |U 

Creutz ratios 20] for the SU(2) and U(l) string ten- 
sions are calculated from Wilson loops up to size 5x5. 
Fig.|3]shows for /3b = 2.3 the behavior of the square roots 
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FIG. 2: Plaquette expectation values on a 12 4 lattice as func- 
tion of A with p a — 1.1 and fib = 2.3 (ordered o and disordered 
d starts). 
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X 

FIG. 3: SU(2) and U(l) string tensions from Creutz ratios 
on a 12 4 lattice as function of A with f3 a = 1.1 and fib = 2.3 
(disordered d and ordered o starts). 



TABLE I: Plaquette expectation values of Fig. [2] 



A 


U(l) 


SU(2) 


Mixed 


0.30 d 


0.716576 (17) 


0.601961 (29) 


0.040142 (08) 


0.40 d 


0.717373 (17) 


0.602630 (27) 


0.054590 (08) 


0.50 d 


0.719702 (21) 


0.604709 (21) 


0.070845 (08) 


0.60 d 


0.721775 (19) 


0.606398 (22) 


0.087592 (10) 


0.70 d 


0.725563 (21) 


0.609602 (29) 


0.107814 (21) 


0.76 d 


0.836426 (08) 


0.723714 (08) 


0.641590 (14) 


0.78 d 


0.837722 (08) 


0.725532 (09) 


0.644591 (14) 


0.80 d 


0.839750 (07) 


0.727838 (09) 


0.650849 (10) 


0.82 d 


0.841460 (34) 


0.729742 (11) 


0.654499 (56) 


0.90 d 


0.850887 (07) 


0.740039 (07) 


0.687311 (11) 


1.00 d 


0.859572 (07) 


0.750205 (08) 


0.712059 (07) 


0.30 o 


0.718566 (19) 


0.603975 (21) 


0.042148 (06) 


0.40 o 


0.719359 (19) 


0.603618 (25) 


0.056604 (09) 


0.50 o 


0.794571 (11) 


0.680028 (10) 


0.461244 (39) 


0.60 o 


0.814685 (11) 


0.700420 (07) 


0.558826 (22) 


0.70 o 


0.829112 (08) 


0.715794 (08) 


0.615873 (14) 


0.80 o 


0.840897 (09) 


0.728681 (09) 


0.656371 (15) 


0.90 o 


0.850874 (06) 


0.740036 (08) 


0.687314 (07) 


1.00 o 


0.859553 (08) 


0.750204 (08) 


0.712054 (09) 



of the string tensions in lattice units a = 1 as function of 
A. Error bars are calculated with respect to 32 jackknifc 
bins, following the scheme of [19]. While the jump in the 
plaquette actions is similar for U(l) and SU(2), this is 
not the case for the string tensions. For SU(2) ^[k, de- 
creases at the transition by a factor 3.5, whereas the drop 
of yfn for U(l) is just 25%. The y/n values are compiled 
in table HIl 

The interpretation of Figj3]is that the U(l) string ten- 
sion signals the deconfined phase on both sides of the 
transition, while the SU(2) string tension is character- 
istic for the confined phase at small A and for a zero- 



TABLE II: String tension values y/n of Fig. [3] 



A 


U(l) d 


SU(2) d 


U(l)o 


SU(2) o 


0.30 


0.0873 (16) 


0.356 (32) 


0.0837 (13) 


0.391 (25) 


0.40 


0.0868 (18) 


0.393 (30) 


0.0855 (16) 


0.426 (22) 


0.50 


0.0835 (16) 


0.418 (18) 


0.07045 (64) 


0.1166 (25) 


0.60 


0.0847 (17) 


0.423 (20) 


0.06694 (45) 


0.1050 (13) 


0.70 


0.0858 (16) 


0.378 (15) 


0.06388 (40) 


0.0982 (13) 


0.76 


0.06117 (35) 


0.0940 (10) 






0.78 


0.11887 (21) 


0.1007 (09) 






0.80 


0.11979 (21) 


0.1027 (10) 


0.06124 (33) 


0.0913 (12) 


0.82 


0.0720 (12) 


0.1051 (10) 






0.90 


0.05839 (32) 


0.0886 (08) 


0.05874 (33) 


0.0869 (07) 


1.00 


0.05714 (28) 


0.0860 (08) 


0.05704 (27) 


0.0847 (09) 



temperature deconfined phase at large A. The latter 
point is supported by comparison with the behavior of 
the U(l) string tension for (3 b = A = at the U(l) decon- 
fining phase transition as shown in Fig. |4] The disconti- 
nuity in the string tension is as in Fig. [3] for SU(2), only 
that no strong metastabilities are observed for U(l). 

Polyakov loop measurements support also the SU(2) 
deconfining phase transition. Fig. [5] shows histograms 
for SU(2) Polyakov loops Pb at A = 0.4 in the confined 
and at A = 0.9 in the deconfined phase. In the con- 
fined phase the values scatter symmetrically about zero, 
whereas in the deconfined phase the Z2 center symmetry 
is broken and the values scatter about a mean of 0.14458 
(45). A very long run would produce a double peak at 
A = 0.9, but our run time was far too short to overcome 
the free energy barrier between the two peaks. Scatter 
plots of the U(l) Polyakov loops at the same couplings 
give in good approximation the ring of Fig. [T]at A = 0.4 
and a more pronounced ring at A = 0.9, both deconfined. 
This interpretation is confirmed by evidence for a mass- 
less photon on both sides and of interest because compact 
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FIG. 4: U(l) string tension from Creutz ratios on a 12 4 lattice 
as function of p a (Pb = A = 0). 
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FIG. 5: Normalized SU(2) Polyakov loop histogram H at 
Pa = 1.1, Pb = 2.3 with A = 0.4 (left) and A = 0.9 (right). 



U(l) LGT allows in contrast to non-compact U(l) LGT 
also for a confined phase. 



C. Spectrum calculations 

Mass spectrum calculations were performed on lattices 
of size N 3 N t , N t > N in the range 4 3 16 to 12 3 64 with 
a statistics of n cq sweeps for reaching equilibrium and 
32 x n cq sweeps with measurements. Data analysis with 
res pec t to 32 bins follows again the jackknife scheme 
of [19| . The n cq values used for different lattice sizes 
are contained in table IIII1 All simulations reported in 
this section are at f3 a = 1.1 and ^ = 2.3. 

First, we estimate the photon mass at A = 0.4 in the 
disordered SU(2) phase and in the ordered SU(2) phase 
at A = 0.9. Following [2l[ this is done via the lattice 
dispersion relation (derived for a free field, for instance, 



TABLE III: Lattice sizes, statistics and photon mass esti- 
mates at A = Ai = 0.4 and A = A2 = 0.9. 



Lattice 




m photon(^l) 


m photon(^2) 


4sin 2 (fci/2) 


4 3 16 


2 11 


-0.252 (20) 


-0.289 (21) 


2 


6 3 24 


2 13 


-0.080 (11) 


-0.0728 (78) 


1 


8 3 32 


2 14 


-0.0138 (18) 


-0.0248 (14) 


0.585786 


10 3 48 


2 14 


-0.0100 (20) 


-0.0123 (11) 


0.381966 


12 3 64 


2 14 


-0.0034 (12) 


-0.0061 (07) 


0.267949 



o 



0.1 



0.01 



0.001 




FIG. 6: Correlation functions data and fits for the photon 
mass estimates. The up-down order of the curves agrees with 
that of the labeling. The upper in a pair of curves is always 
at A = 0.9 and the lower at A = 0.4. 



in[T3) 



in 



photon 



K - 4 sin 2 



(59) 



where E^ is the energy of the U(l) momentum eigenstate 
(fci,0, 0), fci = 2ir/N, of the plaquette trial operator in 
the representation of the cubic group. Energies are 
estimated from correlation functions c(t) of operators by 
the usual two parameter (a\ and E) cosh fits 



z(t) = ax (t 



-Et 



+ e -E(N t -t) 



(60) 



in a range t\ < t < so that the goodness of fit in the 
sense discussed for correlated data in [li| is acceptable. 
Our Ek 1 correlations functions at A = 0.4 and A = 0.9 
are shown in Fig.[B] The corresponding TOp hoton estimates 
are compiled in table Mil The last column of the table 
gives the lowest non-zero lattice momentum squared. 

While in [2l| the photon mass estimates from simula- 
tions on 4 3 16 lattices were within statistical errors con- 
sistent with zero, this is due to increased statistical accu- 
racy no longer true, instead "photon comes ou t negative. 
This holds not just for our action, but also for pure U(l) 
LGT, where we obtained m 2 hoton = -0.2102 (17) on a 
4 3 16 lattice at f3 a = 1.1. Apparently we are dealing with 
an infrared cutoff effect, which disappears with increasing 
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FIG. 7: Glueball and vector boson mass estimates on a 4 16 
lattice (o ordered and d disordered SU(2) starts). 

lattice size and becomes thus consistent with 122j where 
TO photon = was a g am found within statistical errors for 
pure U(l) LGT. So, we conclude that wip hoton — > for 
N — » oo holds also in our cases. 

For an overview of the SU (2) glueball spectrum a 4 3 16 
lattice appears to be large enough, because glueball mass 
values turn out to be high. Compared to the corre- 
lations, the signal from glueball correlation functions is 
very noisy, strongest for the lowest lying zero momen- 
tum A~^ + representation of the cubic group. In Fig. [7J 
we give for the A± + plaquette operator estimates m(l) 
from correlations < t < 1 and m(2) from 1 < t < 2. In 
the disordered phase one has also signals for t = 3 which 
tend to give somewhat lower values than m(2), but are 
so noisy that they overlap within statistical errors with 
m(2). In the ordered phase signals are even weaker, so 
that the correlations at t = 2 include in most cases zero 
within two error bars. We give therefore only the m(l) 
values, which lie considerably higher than in the disor- 
dered phase, with a hysteresis visible at the transition. 
One may conjecture that there are no glueball states in an 
eventual quantum continuum limit of the ordered phase. 
Instead, as shown next, with vanishing SU(2) string ten- 
sion the glueball spectrum appears to break up into mas- 
sive vector bosons. 

Within Higgs model simulations on the lattice trial op- 
erators for the W mass are given by [23J 

W^(x) = -iTt [nW^x)} , (61) 

where W^{x) is the gauge invariant link variable 

W^x) = g\x + fiajV^x) g{x) (62) 

and the SU(2) matrix g(x) collects the angular variables 
of a complex doublet scalar held. This is also a gauge 
invariant operator for our action (l2"TT) with scalar helds. 
Our simulations correspond to this action at very large 
k after fixing the gauge, so that (fUTj) becomes 





A = 0.4 


A = 0.9 


Lat 


m(l) 


m(2) 


m(l) 


m(2) 


mw 


4 3 16 


2.144 (14) 


1.551 (49) 


3.471 (40) 


noise 


0.4706 (65) 


6 3 24 


2.0925 (46) 


1.451 (23) 


3.380 (21) 


2.99 (32) 


0.3385 (32) 


8 3 32 


2.1175 (37) 


1.509 (18) 


3.396 (16) 


2.70 (26) 


0.2901 (13) 


10 3 48 


2.1144 (38) 


1.498 (13) 


3.393 (12) 


2.91 (24) 


0.2729 (09) 


12 3 64 


2.1173 (32) 


1.509 (12) 


3.400 (12) 


2.83 (20) 


0.2677 (09) 




iTr [nV^x)] 



(63) 



1.2 



FIG. 8: Vector boson correlation function on a 4 16 lattice at 
t = 1 (o ordered and d disordered SU(2) starts). 



In the following we calculate vector boson masses from 
correlations of this operator. 

In contrast to the zero momentum correlations of our 
glueball trial operators, those of Vi^ turn out to be beau- 
tifully strong in the ordered phase, while they are zero in 
the disordered phase. For the 4 3 16 lattice the hysteresis 
of the t = 1 correlation is shown in Fig. El At large A val- 
ues (A = 1.0 and A = 1.2) domain walls appear to prevent 
equilibration of disordered starts. Mass estimates in the 
ordered phase from cosh fits (|60[) are found in the lower 
right part of Fig. [7J (left of the transition mw masses are 
infinite) . As these masses rely on long range correlations 
finite size corrections are substantial. For several lattice 
sizes mw estimates at A = 0.9 are collected in table IIV1 
where we also see that there are almost no finite size 
effects for glueball masses. 

The correlation function for the mw fit at A = 0.9 on 
our largest lattice is depicted in Fig. [9] The signal can 
be followed up to more than 20 lattice spacings in the 
N t extension. While significant for small volumes, finite 
size corrections decrease quickly for larger volumes. For 
our N > 8 lattices a fit of the form mw(N) — mw + 
const /N 3, works well and is depicted in Fig. [101 It yields 
the infinite volume extrapolation mw = 0.2567 (10) with 
Q = 0.22 goodness of fit. 
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FIG. 9: Correlation function data and fit for mw mass on a 
f2 3 64 lattice. 
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FIG. 10: Finite size extrapolation of the W mass. 

IV. SUMMARY, OUTLOOK AND 
CONCLUSIONS 

We have introduced a model of SU(2) and U(l) vector 
fields, which can be obtained from the gauge invariant 
interaction (PHI) between a scalar matrix field and the 
vector fields in the London limit (<J> t <I>) 2 — > t (t 2 x 2 
unit matrix), i.e., n — > oo in Eq. (|2ip . On a finite lattice 
this interaction is for sufficiently large k supposed to be 
indistinguishable from the k — > oo limit. MC simulations 
of the quantum field theory on the lattice exhibits a phase 
transition between a deconfined phase with a glueball 
spectrum and a deconfined phase with a massive vector 
boson triplet. A massless photon is found in the spectrum 
of both phases. 

Whether this quantum field theory just lives on the 
lattice or has a quantum continuum limit remains to be 
clarified. The question of its renormalizability requires 
perturbative investigations, which are beyond the scope 
of this paper. Due to the London limit the scalar field 
becomes that of a non-linear a model, which is usually 



non-renormalizable. However, our situation is peculiar, 
because the scalar field can be absorbed into extended 
gauge transformations of the SU(2) and U(l) vector fields 
([T7]). (fT5|) and the theory we are interested in is the one 
without scalar fields as first formulated in [9] . 

In itself the lattice properties are rather remarkable, 
most of all the evidence from Fig. [8] to [10] for a massive 
vector boson triplet in the deconfined, but not in the 
confined, phase. 
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Appendix A: Invariance under Extended Gauge 
Transformations. 

Following .10] we show the invariance of f|2T[) under ex- 
tended gauge transformations. Expanding a calculation 
of [l[ slightly, we find 

9*%K - 9bduB'^ = (Al) 
d^GgaAvG- 1 +i(d u G)G- 1 ] 
~ d^GgbB^G- 1 +i{d tl G)G- 1 ] = 

G(g a d^A„ - gAB^G- 1 (A2) 
+ [{d^G)g a A v - (d u G)gbB l j]G~ 1 
+ GlgaA^G- 1 ) - gtB^G- 1 )} 
+ ip^^G- 1 )-^)^" 1 )] 

Using (d^G-^G + G-^df.G) = d^G^G) = 0, this can 
be transformed to 

QabfX, - g b d v B^ = (A3) 

Gigad^-g^B^G- 1 
+ G { [G~ 1 (d f _ l G),g a A„] - [G^G), g b B^] } G" 1 
- iGiid^G- 1 )^) - {d v G- x )(d ll G)\G-' s - . 

The commutator term transforms as 

ig a g h [B^A' v ]= (A4) 
ig a g h [(GB^G- 1 + (»/ ' g b ){d»G)G- 1 ), 
(GAM- 1 + {i/g a ){d v G)G- 1 )] 
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= ig a9b G[B^A v }G- 1 (A5) 
- GUG-^d^gM - [G-^d^gtB^yG- 1 
+ iGlid^G- 1 )^) - (9„G- 1 )(9 AI G)]G- 1 . 

Combining ([£3]) and (fA"5|) yields (g6|. 



B,(a; + eA) = B„(a:) + ed^B u (x) , (B7) 

are done. These expansions generate more terms (the 
computer does not care) than the analogue expansions of 
(|21[) . but avoid some complications in the identification 
of covariant operators. We write now for i = 1, 2 



Appendix B: Classical continuum limit of the action 
with scalar fields. 



We write the scalar matrix field of the action ([21]) as 
$(x) = 0i(x)0 2 (z) (Bl) 



with cf>i and 02 defined by (|47|) and factor (|2Tj) into the 
form 



S'p, = ^ ReTr (S 1/4 „) Tr (S 2fiv ) + h.c. 
+ rtTr^^-ro) 2 ] . 



(B2) 



As calculation of the classical continuum limit of this ac- 
tion leads to tedious algebra, we rely on the algebraic 
program FORM [24|. For this it is convenient to write 
(|B2[) in a symmetric form with the position x at the cen- 
ter of the plaquette: 



4( 



(B3) 



a 

x — 

a 



x — u — h v— 
F 2 2 



) C/* (* + *|) 0i ( 



a a 
i (x + ^i- + it- 



s' 



(B4) 



.a .a 
2 (z + /z- + 2/- 



^2 ""S J M* - ^ 



'2 I a; - 



2 



, i = 1, 2, to order a 2 . With 
4>i(x + ejl + rjv) = <t>i(x) + ed li (f>i(x) +rjd v (j)i (x) 



We expand the scalar fields 
e and rj of order a: 



+ rjed v d^(j)t(x) + 



(B5) 



The gauge matrices defined by ([7]) have to be expanded 
up to order a 4 in the lattice spacing. For the gauge po- 
tential the substitutions 



S illv ={4>\^) + a 2 S t ^(2) + a 4 S t (4) + . . . . (B8) 

The traces of the order a and a 3 terms of this expansion 
are seen to vanish. The covariant contribution of 

ReTr [S lia ,(2)] Tr [S 2 ^{2)\ (B9) 

to (IB2D comes from the anti-Hcrmitcan terms 



Tr [S 2 ^(2)} = 1 Tr {(£>£ 2 ) f 2 02)* 02 

+ (D* 02 ) 1 02 2 !>* 02 + 2 ) f 02 (D b u 2 ) f 02 

- (l>* 2 ) f 02 02 02 - 4 Dl 02 (l>* 2 ) f 02 



-ft Dl 02 02 02 + 02 Dl 02 (X>° 02 ) 02 

-4 Dl cf> 2 4 D b u ct> 2 -2 4 02 (D b v 2 ) f 02 

-2 05, 02 (l>* Dl 2 ) f 02 + 2 2 02 (D» 2 ) 1 £>£ 02 
+ 2 0^02 02^^02} (BIO) 

and Tr [S'i /J j/(2)], which is obtained by interchanging A a 
with B a {D% with D h a ) and then the subscripts /i and v. 

In the London limit Tr [Si^(2)] Tr [5 2/1I /(2)] gives the 
contribution ~ Tr (d^A^d^B^) to the action (j2"T)) . The 
terms - Tr (ft,j4 M ft,^) and ~ Tr (9^ fyB,,) of (gft 
come from the a 4 (0| 0i) 2 5 2/J! /(4) and a 4 5i^(4) (0| 2 ) 2 
contributions to (IB2|) , whose calculation we found too te- 
dious to pursue. For small oscillations about the expec- 
tation value of the scalar fields 

5i^(4) = ^(2) 2 and S 2 „„(4) = S 2 ^{2) 2 (Bll) 

are good approximations, illustrating the typical gauge 
invariant terms encountered. 



A M (x + T)i>) = A M (x) +rjd u A fl (x) , 



(B6) 
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